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AN IRREVERSIBLE LOCAL MARKOV CHAIN THAT PRESERVES 
THE SIX VERTEX MODEL ON A TORUS 

ALEXEI BORODIN AND ALEXEY BUFETOV 


Abstract. We construct an irreversible local Markov dynamics on configurations of 
up-right paths on a discrete two-dimensional torus, that preserves the Gibbs measures 
for the six vertex model. An additional feature of the dynamics is a conjecturally 
nontrivial drift of the height function. 


1. Introduction 

Random growth models is a rapidly developing subject that focuses on studying 
large-time behavior of randomly growing interfaces in (d+l) dimensions. The growth 
mechanism is usually assumed to be local, in the sense that distant parts of the interface 
evolve (almost) independently, to conform with the belief that most inter-molecular 
interaction mechanisms in nature are local. 

Not surprisingly, the best understood case is d = 1. Many results are available includ¬ 
ing the hydrodynamic (law of large numbers) behavior for certain classes of systems, 
see e.g. the book of Kipnis-Landim Q and reference therein, as well as fluctuations 
for certain integrable systems in the KPZ universality class, see e.g. Ferrari-Spohn 
[To], Corwin [8], and references therein. 

One key fact that is very useful in the (l-l-l)-dimensional situation is that the steady 
states of the growth model often turn out to be described by product-measures. The 
lack of such a simple structure in higher dimensions is one obstacle for studying the 
d > 1 case. Furthermore, it is a challenge to find any irreversible (or driven) Markov 
evolution of an interface in (d-i- l)-dimensions with d > 1 for which one could describe 
the steady states. See Gates-Westcott [HI for a notable exception. 

Recent developments in studies of integrable probabilistic models brought up new 
examples. In i, Borodin-Ferrari constructed a local random growth model of 
stepped surfaces in (2 -i- l)d for which the asymptotic growth velocity (as well as 
local correlations and global Gaussian fluctuations) were computed explicitly, see also 
Borodin-Bufetov-Olshanski [3] for a generalization that includes similar results for 
various initial conditions. 

Very recently, Toninelli [TS] proved that the dynamics of Borodin-Ferrari, as well as 
its generalization that acts on dimers on the square lattice, can be correctly defined 
on translation-invariant random surfaces over the whole plane. He also proved that 
members of a particular two-dimensional family of such random surfaces are steady 
states for the resulting dynamics. In fact, this two-dimensional family exactly coincides 
with the two-dimensional family of translation-invariant Gibbs measures for dimers on 
either hexagonal or square lattice, cf. Kenyon-Okounkov-Shelfield [12], Sheffield [IS]. 
While checking similar statements on the torus in this case is fairly straightforward, 
extending to the whole plane is very nontrivial and requires substantial efforts. 

^We omit here the case of heat bath or Glauber-type dynamics that are typically reversible, have 
asymptotically zero drift, and belong to different universality classes than the irreversible examples. 
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The correlation functions of the measures from this two-dimensional family are deter- 
minantal. These measures can be thought of as originating from free-fermion (dimer) 
models. In the last few years certain progress has been achieved in understanding 
random growth models in (l+l)d that are not free fermionic (but integrable), see e.g. 
[7] and references therein. It is enticing to try to expand this progress into the domain 
of (2-|-l)-dimensional models. 

One step in this direction was taken by Corwin-Toninelli [9] . They used a triangular 
array dynamics from Borodin-Corwin [Ij (which generalized [5] and [2] to a non-free 
fermion case) to define an irreversible Markov dynamics on particles living on a two- 
dimensional torus. They found that Gibbs measures with respect to a particular 
explicit g-weighting are preserved by the dynamics. 

In this work we do something similar, but the extension of 0.0 is different, and 
it is not a part of the Macdonald processes framework of [Ij . 

We construct an explicit irreversible Markov dynamics on up-right path configura¬ 
tions on a discrete two-dimensional torus. Such path configurations can be viewed as 
level lines of a height function whose plot is the evolving interface. Our main result 
states that the Gibbs measures of the celebrated six vertex model are invariant with 
respect to our dynamics. 

The literature on the six vertex model is vast, and we will not even try to survey it. 
An interested reader could consult e.g. Palamarchuk-Reshetikhin [H] and references 
therein. Its Gibbs measures on a torus are parametrized by two nonnegative integers 
(in addition to the weights of the six vertex model itself). 

There remains an open and very interesting question to see if the dynamics defined 
in this paper can be extended to the whole plane while preserving the two-dimensional 
family of the six vertex translation invariant Gibbs measures. Many things are known 
about such Gibbs measures but most of them are conjectural, although widely accepted 
in the physics literature. One might hope that having an explicit dynamics that 
preserves them could shed more light on their properties. 

In the text below we give a verification proof of our main result. A brief explanation 
of its derivation is given in the Appendix; it is based on properties of a family of 
rational symmetric functions introduced in [T]. Guriously, while the positivity of the 
measures participating in the derivation requires that the parameters of the six vertex 
model he in the so-called ferro-electric region, the final result knows nothing about 
this restriction — the Markov dynamics is positive and preserves the Gibbs measures 
for arbitrary nonnegative values of the six vertex weights of the model. 

Acknowledgments. A. Borodin was partially supported by the NSF grant DMS- 
1056390. A. Bufetov was partially supported by Dynasty foundation, the RFBR grant 
13-01-12449, and by the Government of the Russian Federation within the framework 
of the implementation of the 5-100 Programme Roadmap of the National Research 
University Higher School of Economics. 

2. Preliminaries 

Let Th be the square grid on a two-dimensional torus of size M x N. It can be rep¬ 
resented by the (conventional) grid {0,1,..., M} x (0,1,..., N} with the identification 
of vertices (0,f) = {M,i) and (qO) = z = 0,1,.... We say that two vertices are 

adjacent if they connected by an edge in Th; we say that the vertices are h-adjacent 
if this edge is horizontal (the first coordinate changes) and v-adjacent if this edge is 
vertical. We will use the circular order in the horizontal and vertical directions that is 
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Figure 1. We denote these 6 types of vertices by symbols ©, +, r, j, 
I, —, respectively. 



inherited from the sets and All relevant operations 

and notions (like “a + 1” or “the top vertex from v-adjacent vertices”, etc) are taken 
in accordance with these circular orders. 

A state of a 6-vertex model on Th is a subset of edges from Th such that each vertex 
has one of 6 allowed types. We will think of a horizontal edge from this subset as an 
arrow directed to the right and of a vertical edge as an arrow directed to the top. The 
6 allowed types of vertices and our notation for them are shown in Figure [!} and an 
example of a state is shown in Figure To each of these types of vertices we assign a 
real positive weight] we denote these numbers byr<;(0), w{+), w{r), tc(|), and 

w {—). The weight of a state is defined as the product of weights of all vertices from 
Th. 

For a state s, let Ns{X) denote the number of vertices of the type X in s (here 
X can be one of 6 allowed types of vertices). For a vertex A from Th we denote by 
types(A) the type of this vertex in the state s. 

It is clear that the number of horizontal arrows in any state equals Mki , for some 
fci € N u {0}, and the number of vertical edges equals Nk 2 , for some ^2 e N u {0}. Let 
Ski,k 2 be the set of all states with fixed numbers ki,k 2 . The Gibbs measure on 
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Figure 3. Jumps to the right happen when there are two vertices with 
such a structure of arrows; arrows on edges not in the picture can be 
arbitrary. The four different cases are shown in Figures]^ — 


is defined by 


Prob(s) 


w{s) 


^ ^ ^ki,k2 ■ 


If all vertex weights are positive, then this is a probability measure on Sk^^k 2 - Though 
all 6 weights of vertices can be arbitrary positive numbers, this measure depends on 
only two free parameters; see the end of Section]^ for more details. 


3. Description of the dynamics 


The goal of this section is to describe our dynamics V. This will be a continuous 
time Markov chain on the f ini te set We will consider the case l<A:i<A^-l,l< 

k 2 < M -1 only. In the remaining cases fci = 0, or A :2 = 0, M the set of states becomes 
rather poor, and our dynamics degenerates to a one-dimensional dynamics (ASEP) on 
it. 

For the description it is enough to define the matrix elements p(si, S 2 ), Si, S 2 ^ Sk^^k 2 
of the generator of our Markov dynamics. Since ki,k 2 are fixed, we will omit the 
dependence on them in the sequel. 

Let us proceed with a formal description. A jump to the right can occur when 
there is a configuration of arrows depicted in the left part of Figure]^ (here and in the 
sequel we assume that arrows on the edges that are not present in the figures can be 
arbitrary). We want to single out all types that vertices A and B from Figurecan 
have. Let us call the pairs of h-adjacent vertices depicted in the left panels of Figures 
S-EI r-admissible pairs. 

In a similar vein, we call the pairs of h-adjacent vertices that are depicted on the 
left panels of Figures I-admissible pairs. 

Let us define the possible transform of a state s e S; there are two types of such 
transforms — for jumps to the right and for jumps to the left. 

Let us start with describing the possible transforms that correspond to the right 
jumps. Let A and B be an h-adjacent pair of vertices from Figures i-0 Let us 
denote by C the vertex at the same vertical line as A such that all vertices between 
A and C have type | and C has type r or +. Note that such a vertex always exist 
since fci > 1 by our assumption, i.e. there should be a horisontal arrow in s that starts 
on the vertical line which contains A. Let D be the vertex that is h-adjacent to C 


and is in the same vertical line as B (see the left panel of Figure 12 for the pictorial 
explanation of these definitions). 
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Figure 4. Admissible pair; jump to the right; rate = 




w(+)w((Z)) ‘ 





n 




Figure 5. Admissible pair: jump to the right; rate = , , ^ • 


The possible transform corresponding to the jump to the right changes the edges in 
the rectangle ABCD as shown in Figure the state of all other edges remains the 
same. 

Let us now describe the possible transforms that correspond to the jumps to the 
Let A and B be an h-adjacent pair of vertices from Figures — pT 


left. 


Denote 


by C the vertex at the same vertical line as A such that all vertices between A and 
C have the type — and C has the types © or J. Note that such a vertex always 
exist since fci < A” - 1 by our assumption, i.e. there should be a horisontal edge free of 
arrows of s that ends on the vertical line which contains A. Let D be the vertex that 
is h-adjacent to C and is in the same vertical line as B (see the left panel of Figure [l3| 
for the pictorial explanation of these dehnitions). 

The possible transform corresponding to the jump to the left changes the edges in 
the rectangle ABCD as shown in Figure [T^ the state of all other edges remains the 
same. 

Now we are in a position to define our Markov dynamics V. The transition rate 
of our dynamics p{si S 2 ), -51,52 ^ S, is positive if and only if one can chose an 
admissible pair of vertices in si such that the corresponding possible transform turns 
Si into S 2 . The value of p(si S 2 ) depends only on the type of admissible pairs and 
is given in Figures]^ — [^for each of them. 

Let us comment on our exact formulas for the transition rates. 

First, note that the r-admissible and /-admissible pairs are in the natural bijection: 
The jump to the right in Figure |^ ( [^ [^ [^ is exactly the reverse to the jump to the 
left on Figure]^ (|^ respectively). The jump rates for corresponding pairs are 
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Figure 6. Admissible pair; jump to the right; rate = ) 




n A 







Figure 7. Admissible pair; jump to the right; rate 
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Figure 8. Admissible pair; jump to the left; rate = 
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Figure 9. Admissible pair; jump to the left; rate 
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Figure 10. Admissible pair; jump to the left; rate 
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Figure 11. Admissible pair; jump to the left; rate 
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Figure 12. Jump to the right results in such a change of vertices. The 
vertices between A and C can have type | only; however, there can be 
arbitrarily many of them. Not depicted arrows can be arbitrary and 
remain the same. 


also inverse to each other. Second, the four different jump rates for the jumps to the 
right are also split into two pairs which are inverses of each other. Note also that the 
values of jump rates depend only on local properties of a state. However, the jumps 
themselves can propagate upwards arbitrarily far. 

The statements of the next two paragraphs will not be used in the sequel, but it is 
natural to address them here. 
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Figure 13. Jump to the left results in such a change of vertices. The 
vertices between A and C can have type | only; however, there can be 
arbitrarily many of them. Not depicted arrows can be arbitrary and 
remain the same. 


There is another bijection between jumps to the right and to the left. For a state s 
one can define a dual state which is defined as the set of all edges except for the edges 
from s. Then the jumps to the left (Figures]^ 10, are in correspondence with 
the jumps of the dual state to the right ( Figures Hlli respectively). The jump 
rates for corresponding pairs coincide. 

The jump rates for any admissible pair can be computed with the use of dual 
states. The recipe is: Multiply the weights of the pair of vertices and their weights 
in the dual state (four factors in total) after the transform and divide by the same 
product computed before the transform; then take the square root. For example, 
for an admissible pair from Figure |4 this computation looks as follows. After the 
transform (the right panel of Figure ^ we need to take w( —)r’(-|) — the weight 
of these vertices — and then multiply it to w{\)w{r) — the weight of these vertices 
in the dual state. Analogously, before the transform (the left panel of Figure we 
obtain w{-i)w((D) x w{r)w{+). Dividing these quantities and taking the square root 
we obtain the jump rate for this admissible pair from the caption of Figure]^ 

We have defined all jump rates for our dynamics; thus, our description of the dy¬ 
namics T> is complete now. Let us add more comments. 

We operate with 6 weights for vertices without any constraints; however, different 
weights can give the same Gibbs measure. We have 4 equalities 

1) Ns{+) + Ns{-\) + Ns {—) = const — because all configurations must have the same 
number of horizontal arrows. 

2) Ns{+) + Ns{-i) + iVs(l) = const — because all configurations must have the same 
number of vertical arrows. 

3) Ns{Q) + Ns{+) + Ns(-i) + Ns(r ) + Ns (—) + NgQ) = const — because the number 
of all vertices is the same for all configurations. 

4) Ns(-i) - Ns{r ) = 0 — see the proof of Lemma below. 

This gives rise to 4 possible types of transforms of weights that produce the same 
Gibbs measure. They are as follows: 

1) Multiply R’(+), 'R'(j), and w {—) by Ci. 

2) Multiply w{+), and t(j(|) by C 2 . 
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Figure 14. An example of a state and its flip on a torus grid. 

3) Multiply the weights of all vertices by C 3 . 

4) Multiply w{-i) by C 4 and multiply w(r) by 

One easily checks that our jump rates are invariant under these transforms; it means 
that for the same Gibbs measure we have the same dynamics. 

We also note that since we have 6 weights for the types of vertices and 4 possible 
transforms of them, then in fact we have only 2 free parameters of the model coming 
from weights; we also have two integer parameters ki and /c 2 . We have no further 
restrictions: If the weights of vertices are positive, then all our jump rates are positive 
as well. 


4. Invariance of the Gibbs measure 

In this section we establish the main theorem of the present paper, Theorem]^ below. 
We start with some necessary notions and lemmas. 

For any state s e we define a state s which we call the flip of s. By definition, this 
is the state which is obtained from s by changing all arrows into the opposite direction 
and by using the opposite circular order in both horizontal and vertical directions. In 
other words, we rotate the grid into 180 degrees and change the direction of all arrows. 
See an example in Figure 

Lemma 1. For any s € S we have w{s) = w{s). 

Proof. Note that vertices of types +, 0, |, and — are invariant under the flip, while 
r turns into J and vice versa. Therefore, it is enough to prove that Ns{r) = Ns{-\) 
for any s e S. But this is clearly visible from the up-right paths structure of s: Along 
each up-right path the number of right turns equals the number of left turns, and each 
vertex + contains one left turn and one right turn inside. □ 

Lemma 2. Let si,S 2 c S. If p(si ->■ S 2 ) > 0, then p(s 2 si) > 0 as well. Moreover, we 
have 

(1) t«(si)p(si ^ S 2 ) ='»^(S2)P(S2 ^ Si). 

Proof. Let us consider the case when si transforms into S 2 by a right jump. 

From the construction, the two top vertices — G and D — in which the transforma¬ 
tion happens can look in 4 possible ways (see Figure [T^ and Figure [l5|. It is readily 
visible that after the flip they exactly coincide with admissible pairs of the jumps to 
the right from S 2 to Si (left panels of Figures]^ respectively). This proves the 

first statement of the lemma. 

















































10 


ALEXEI BORODIN AND ALEXEY BUFETOV 










Figure 15. The 4 possible ways how the non-depicted arrows in C and 
D from Figure [T^ can look like. 


Let us prove now the second claim. Lemma [I] asserts that w(s 2 ) = w{s 2 ). Note that 
the set of types of vertices between A and C, and between B and D remains the same 
during the transform (see Figure [I^. Thus, we have 

. ^ . w(type3^(^))w(type3^(.B))u;(type,^((7))w(type3^(T>)) 

W S2 W 


since we need to control the change of types in A,B,C,D only. Substituting these 
equalities into Q we see that it remains to prove 

p{sx ^ 52 ) ^ '»^(type.2(^))^(^yp^«2(-^))^(^yp^.2(<^))'»^(^yp^«2(-P)) 

p(s 2 ^ si) w(type 3 ^ (24))M;(type,^ (S))M;(type 3 ^ (C'))w(type,^ (D))' 


This equality can be directly checked from the definitions of the jump rates. We give 
this verification in the form of a table below. 


/ p(si^S2) 
p(s2->-Sl) 

Figure]^ 

Figure]^ 

Figure]^ 

1 

«;(—)m)(|) 
w;( j)i(;(r) 

Figure 

w;(j)K;(r) 
w( — )w(\) 

1 

Figure]^ 

w(+)w{(D) 

w{+)w((D)w{ — )w{\) 

w(A)w(r ) 

w(A)‘^wlr )2 

^Figure 

w(+)w((D) 
«;(—)?x)(|) 

w(+)w((D) 

w{A)w{r) 


Figure]^ 

w(+)w(Q) 
w(A)‘^w(r)^ 
w(+)w((D)w( — )w(\) 

1 

k;( j)K;(r) 

?X)(—)«;(!) 


Figure 

«’(—)«’(!) 

w(+)w((Z)) 

w(A)w(r) 

w(+)w() 

«’(—)«’(!) 

k;( j)K;(r) 


1 / 


The rows of this table correspond to different types of vertices in A and B. The 
columns correspond to different types of vertices in C and D after the flip. For the 
entries we simply calculate the quotient of jump rates (which are defined in the caption 
to the corresponding figures). On the other hand, it can be directly checked that the 
numerator of the obtained result corresponds to the weights of vertices which appear 
after the transform, while the denominator correspond to the vertices which disappear 
after the transform. 

For example, if the vertices A and B are as in Figure and the vertices C and D 
after the flip are as in Figure]^ then in si (before the jump) the vertices A, B, C, and 
D have types J, r, r, —, respectively, and in S2 (after the jump) they have types —, 
+, 0, r, respectively. Therefore, the quotient of the product of all weights is equal to 
w(j)w(r) ’ which equals the corresponding quotient of jump rates. 

The verification in other cases is analogous. 

The case of jumps to the left can be proved in a similar way. However, it also follows 
from the case of jumps to the right. Indeed, in the case of jumps to the left all changes 
of types of vertices in A, B, C, D are inverse to the changes that happen with the 
corresponding jump to the right (recall that the r-admissible pairs and /-admissible 
pairs are in the natural bijection). On the other hand, the jump rates for corresponding 
jumps are inverse to each other as well. Therefore, we have the desired equality in the 
case of jumps to the left as well. 


□ 
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Figure 16. An example of configuration in Th;t and blocks inside it. 


For two types of allowed vertices A, Y and a state s, let us denote by Ns(X,Y) the 
number of /i-adjacent pairs such that the left vertex from this pair has type X and the 
right vertex from this pair has type Y. We are interested in these numbers because 
the jump rates of our dynamics are determined by pairs of h-adjacent vertices. 

Lemma 3. For any state s e S, we have 

(2) W( J, O) + w(—, J) = w(0, r) + W(r, —). 

Proof. All quantities in Q depend on two h-adjacent vertices. Recall that our model 
lives on the grid (0, 1) • • •) Af - 1} x {0, !> • • •) “ !}• Let us consider the part of our 

grid Thfc := (fc, A: + 1} x {0, !> • • • > “ 1}) for some 0 < k < M. Let us prove that we have 

the desired equality of pairs of h-adjacent vertices inside Th^. 

We say that i € (0,1) • • •) Ai - 1} is a good number if two vertical edges of the form 
{k,i) (k,i + 1) and {k + l,i) (k + l,i + 1) do not contain arrows from s. We say 
that a block is a collection of good numbers {o, a + 1, o + 2,. .., 6} such that o - 1 and 
6+1 are not good numbers. See Figurefor an example. 

If all numbers from 0 to A -1 are good, then we obviously have the desired equality. 
Otherwise, note that for each block {a, a + 1, a + 2,..., 6} the vertices {k, 6 + 1) and 
(A: + 1,6 + 1) must be of the form (j, 0) or (—, j). On the other hand, the vertices 
{k,a) and {k + l,a) must be of the form (0,r) or (r,—). Therefore, the number 
of blocks equals the number of configurations both in the left-hand side and in the 









12 


ALEXEI BORODIN AND ALEXEY BUFETOV 


right-hand side of Q. Thus, we have this equality inside each Th^. Summing over all 
k, we arrive at the statement of the lemma. □ 

After the flip of a state, a pair of h-adjacent vertices transforms into another pair of 
h-adjacent vertices. The next lemma compares the numbers of a certain pairs of types 
of /i-adjacent vertices in a state and in its flip. 

Lemma 4. Let 

A:=iV,(j,©)-iV,(©,r). 

Then we have 

(3) Ns{\,r) - NsU I) = A, Ns{Q, |) - iV^d, ©) = A, 

and 

(4) iV,(r,—,j) = A, 7V,(+,j)-iV,(r,+)=A, iV,(—,+)-iV,(+, —) = A. 

Proof. In the proof of Lemma above we showed that = Ns(r). We have 

iV,( J) = N,U ©) + N,{J, I) + N,{J, r), 

A^,(r) = A^3(©,r) +A^3(|,r) +A^3(j,r). 

Subtracting the second equality from the flst one, we obtain A^s(|,r) - A^s(j, |) = A. 
Similarly, 

W(|) = iV.(O,|) + iV.0,|) + iV.(|,|), 

W(|) = iV.(|,©) + iV.(|,r) + iV.(|,|). 

Subtracting the second equality from the first one we obtain the remaining equality 
from (|^. Applying already proved equalities Q to the dual to s state, we obtain that 

iV,(r, — ) - Ar,(— , j) = S, 7V,(+, j) - 7V,( j, +) = B, 7V,(—, +) - iV,(+, —) = S, 

for some B eZ. Finally, Lemma 1^ gives A = B. □ 

Now we are ready to prove our main result. 

Theorem 5. The Markov dynamies V preserves the Gibbs measure for the 6-vertex 
model on a torus. 

Proof. Let be the set of all pairs of states Si, S 2 e such that p(si ^ S 2 ) is strictly 
positive. It is enough to prove that for each s € S we have 

(5) ^ w(s)p(s^si)= ^ w(s 2 )p(s 2 ^ s). 

5i:(s,si)e7^+ S2-(s2,s)eV'^ 

Using Lemma we obtain 

E W(S2)P{S2 ^ S) = ^ w(s)p{s^S2). 

S2-(s2:S)eV'^ S2-{S,S2)^V'^ 

Taking into account Lemma and the equation above, we reduce (|^ to 

(6) X! Pis^si)= X! P(s^S2). 

Si:(s,Si) 67^+ S2:(5,S2)eP+ 

Recall that the jump rates are determined in an explicit and rather simple way, and 
note that Lemma provides us with the required information for comparing the left- 
hand and the right-hand sides of ([^. Indeed, all quantities in Lemma 4 compare the 
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number of the same admissible pairs in a state and its flip. Taking into account all 
admissible pairs, we get 


(7) X! Pis^si)- Y. P(s^S2) 

si:(s,si)€P+ S2-(s,S2)^V'^ 

^ w{—)w{\) ^ w{-i)w(r) 


w(+)w(Q) 


+ 


-A 


)w{Q) ^^w(+)w{Q)w( —)tc(|) \ 

w{- 



w{-\-)w((D) ^ ^ w(-i)w{r) 

■0 

w(—)w(|) 

w( )w(|) ^w(-\-)w{Q)w( —)tc() 

w(+)«;(0)^ 


The first parentheses come from the jumps to the right; the second parentheses come 
from the jumps to the left. Indeed, the first term in the first parentheses comes from 
the comparison of /i-adjacent vertices (-1,0) (Figure]^ in a state and in its flip, the 
second term in the first parentheses comes from (|, 0) (Figure [^, and the third term 
from the first parentheses comes from (|, j) (Figurej^. The pair (J, r) (Figure]^ does 
not contribute because it is invariant under the flip. 

Analogously, the first term in the second parentheses comes from the comparison 
of h-adjacent vertices (—, j) (Figure]^ in a state and in its flip, the second term in 
the second parentheses comes from (—,+) (Figure [To)), and the third term from the 
first parentheses comes from (r,+) (Figure [TT|. The pair (r, j) (Figure]^ does not 
contribute because it is invariant under the flip. 

It is clear that the total sum in ([^ equals 0, which concludes the proof of the 
theorem. □ 


Remark 6. The computations above shows that in the general case we need both 
jumps to the right and to the left in order to obtain the invariance of the Gibbs 
measure. However, at the free-fermion point (i.e., in the case w{-i)w{r)-w{+)w{Q) - 
w{ —)'ii'(|) = 0) the dynamics can be decomposed into two parts; the part of the 
dynamics that corresponds to jumps only to the right still preserves the Gibbs measure 
(and the part of the dynamics that only jumps to the left as well). This is clearly visible 
from ([^ — both parentheses vanish. 


5. Appendix: The origin of dynamics 


In this section we briefly describe the origin of the dynamics V. 

Let u and q be two reals, and consider the following weights of the 6 types of vertices: 


(8) tc(0) = 1, 'w{+) = 


u 


\/Q 


1 - 

w(|) = 


w(—) = 


u ■ 


1 


1 _ ’ 


1 - 
w{-i) = 


Ai 

1-g 

1 _ ^ ’ 
A 


w(r) = 


(1-9 

1 - ^ 

A 


These weights coincide with weights from pQ Definition 2.1] with s = q~^R. 

Let A = (A(l), A(2),..., X{N)), A(l) > A(2) > • • • > X{N), be a A^-tuple of decreasing 
nonnegative integers. Let Sign^ be the set of all such A^-tuples (A’s are called signa¬ 
tures). We say that A^v ^ Sign^ and A^v+i ^ Sign^^;^ interlace and write A^v < Aat+i 
if 


Ajv+i(l) ^ Ajv(l) > Ajv+i(2) > • • • > X]\f{N) > Ajv+i(A^ + 1). 
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ALEXEI BORODIN AND ALEXEY BUFETOV 


Define a triangular array of particles as the collection Ai < A 2 < • • • < Ajv (it is usually 
referred to as a Gelfand-Tsetlin scheme due to its relation to the branching of the 
irreducible representations of the unitary groups). 

For A € Sign^ let Fx and Gx be the rational symmetric functions introduced in [T]. 
They obey an analog of the Cauchy identity m Equation (4.6)] and its skew versions). 
This allows us to apply a general construction of dynamics on triangular arrays which 
samples from a certain probability measure. The idea was introduced in [5], [2] for 
the Schur symmetric polynomials and later applied in [H Section 2.3] to Macdonald 
polynomials. The exposition in [H Section 2.3] can be applied to the case of F- and 
G- functions directly. 

One of the main properties of such a dynamics is that it samples from a probability 
measure on a triangular array. Let K he a, positive integer which will play a role of 
time in our dynamics, and let x be a real fixed parameter. It is possible to sample a 
probability measure such that for each N the projection to Sign^ of this measure has 
the form 

const ■ Fx(x, X,... ,x)Gx(u, u,... ,u) 

N K 

(in fact, the joint probabilities for different N also have a prescribed form which is 
similar to the ascending Macdonald process, see [H Section 1.1]). Moreover, this 
dynamics allows to make a step K ^ K + 1 (that is, starting with the probability 
measure corresponding to K obtain the probability measure corresponding to K + 1) 
by rules which involve only local interaction of particles of the triangular array. We 
denote by A^ the (random) signature from Sign^ at time K. 

Due to the combinatorial formula for F- functions (see [5], Definition 3.1]) this dy¬ 
namics has a rather simple form. Let us further assume that u = q~^G + s, for 0 < £ « 1. 
Then the vertices | and — have small weights proportional to e (see Q). This im¬ 
plies that with large probability the step K ^ K + 1 oi our dynamics on each level 
looks like A|^'^^(f) = A^(f) + 1, i = 1,2 ,..., N, and the events A|^'^^(f) = X^{i) and 
A^'^^(z) = \%{i) + 2 happen with probability proportional to s, while all other events 
can be considered as impossible (these two ^-probability events eventually give rise 
to the left and right jumps in dynamics V). We are interested in the limit of this 
dynamics for N = 0. If we consider this dynamics “in the bulk”, then our 

formulas for probability measures (heuristically) imply the preservation of the Gibbs 
measure on the set of paths. The dynamics on the triangular array is governed by local 
interactions only. Therefore, it is reasonable to assume that “boundary conditions” 
do not affect the dynamics in the bulk, and that after the limit £ ^ 0 the dynamics in 
the bulk will preserve the Gibbs measure on a six-vertex model on torus. This paper 
describes the limit dynamics and proves that it has such a property. 
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